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We  present  a  zero  dimensional  stationary  model  which  precisely  predicts  the  current-voltage- 
characteristics  of  anode  supported  SOFC  single  cells  over  a  wide  range  of  operating  conditions.  The 
different  kinds  of  electrode  polarization  resistances  are  separated  from  experimental  impedance  data 
by  means  of  a  detailed  equivalent  circuit  model  developed  specifically  for  the  analyzed  cell  type.  The 
activation  losses  are  modeled  by  the  Butler-Volmer  equation,  whereas  the  loss  contributions  from  gas 
diffusion  polarizations  are  calculated  from  Fields  law.  The  partial  pressure  and  temperature  dependency 
of  the  cathodic  and  anodic  exchange  current  density  could  be  determined  by  a  fit  of  semi  empirical  power 
law  model  equations.  The  exponents  c  and  d  for  the  CO  and  C02  partial  pressure  dependency  of  the  anodic 
exchange  current  density  are  determined  independently  of  each  other.  This  paper  presents  the  modeling 
results  for  a  wide  range  of  operation  parameters  as  well  as  their  experimental  verification. 

©  2010  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Interest  in  solid  oxide  fuel  cells  (SOFC),  which  combine  poten¬ 
tially  efficient  energy  production  with  flexible  fuel  composition, 
has  significantly  increased  in  the  past  few  years.  Amongst  the 
most  important  development  goals  is,  besides  long-term  stabil¬ 
ity,  the  single-cell  performance,  which  is  limited  by  three  loss 
mechanisms:  (i)  ohmic  losses,  (ii)  polarization  losses  through  gas 
diffusion,  and  (iii)  losses  by  activation  polarization. 

In  this  study  the  different  loss  mechanisms  and  their  contri¬ 
butions  are  identified  and  subsequently  separated  by  means  of 
impedance  spectroscopy  in  combination  with  a  fit  of  the  mea¬ 
sured  spectra  to  a  detailed  equivalent  circuit  model  developed 
specifically  for  the  cell  type  investigated  and  presented  in  Ref. 
[1].  The  parameters  determined  in  this  manner  are  used  to  pre¬ 
dict  the  current-voltage  ( l-U )  characteristics  of  the  analyzed  anode 
supported  SOFC  single  cells  (ASCs)  in  a  broad  range  of  operating 
conditions  [2]  (Fig.  1). 

2.  Experimental 

The  ASCs  investigated  are  based  on  50  mm  x  50  mm  anode 
substrates  (Ni/8YSZ,  thickness  1  mm),  an  anode  functional 
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layer  (Ni/8YSZ,  thickness  approximately  7  pim),  a  thin  film 
electrolyte  (8YSZ,  thickness  approximately  lOpim)  and  a 
Ce0.8Gdo.202-5  (CGO)  interlayer  (thickness  approximately  10  pm). 
The  La0.58Sro.4Coo.2Feo.803_5  cathodes  (lcmxlcm,  thickness 
approximately  45  pm)  were  deposited  by  screen  printing  and 
subsequently  sintered.  Detailed  information  on  preparation  and 
properties  of  these  cells  is  given  in  Refs.  [2-5]. 

The  cells  were  characterized  by  electrochemical  impedance 
spectroscopy  (EIS)  in  a  wide  range  of  operating  conditions  (vari¬ 
ous  temperatures  and  partial  pressures  of  carbon  monoxide  and 
carbon  dioxide  (pCO,  pC02)).  The  EIS  measurements  were  carried 
out  with  a  Solartron  1260  FRA  in  the  frequency  range  from  30  mFIz 
to  1  MHz. 

l-U  measurements  were  carried  out  at  temperatures  between 
650  °C  and  800  °C  with  current  densities  up  to  2  A  cm-2  and  cell 
voltages  not  below  a  limit  of  0.6  V.  This  operation  range  guaranties 
a  stable  cell  performance  during  all  experiments. 

The  single  cells  were  investigated  under  ambient  air  pressure 
at  various  gas  compositions  of  CO,  C02,  and  N2  at  the  anode  side. 
Pure  air  was  used  on  the  cathode  side.  The  overall  gas  flow  rate  was 
constantly  kept  at  250  ml  min-1  for  the  anode  as  well  as  the  cath¬ 
ode.  These  fairly  high  flow  rates  were  chosen  to  assure  constant 
gas  compositions,  i.e.  a  “low”  fuel  utilization  within  the  entire  cell 
area.  For  cell  voltages  not  lower  than  0.6  V  a  maximum  fuel  uti¬ 
lization  of  about  9%  can  be  calculated.  For  the  analyzed  operating 
conditions  the  latter  corresponds  to  a  maximum  voltage  drop  due 
to  fuel  utilization  not  higher  than  10  mV  (calculated  by  the  Nernst 
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Fig.  1.  Equivalent  circuit  model  used  for  the  CNLS  fit  of  impedance  data  [1  ]. 


equation).  For  this  reason  in  the  present  work  the  voltage  drop  due 
to  gas  conversion  has  been  neglected. 

3.  Evaluation  of  impedance  spectra 

The  measured  impedance  curves  were  evaluated  using  the  com¬ 
mercially  available  complex  nonlinear  least  squares  (CNLS)  fit 
software  Z-View  [6].  The  equivalent-circuit  model  used  for  this 
procedure  (Fig.  1)  was  developed  by  a  pre-identification  of  the 
impedance  response  by  calculating  and  analyzing  the  correspond¬ 
ing  distribution  of  relaxation  times  (DRT)  and  has  been  presented  in 
Refs.  [1,2, 7-9].  For  a  more  detailed  description  of  the  DRT  method 
and  its  application  we  refer  the  reader  to  Ref.  [10]. 

The  applied  equivalent  circuit  (Fig.  1)  consists  of  5  serial 
impedance  elements:  (i)  an  ohmic  resistor  (R0-  ohmic  losses),  (ii) 
two  serial  RQ  elements  (R2A  and  R3A:  gas  diffusion  and  activa¬ 
tion  polarization  in  the  anode  functional  layer  [1,7,9,11]),  (iii)  a 
Gerischer  element  [12]  (R2C:  activation  polarization  of  the  cath¬ 
ode),  (iv)  a  generalized  finite  length  Warburg  Element  (G-FWS) 
(R-ia-  gas  diffusion  polarization  in  the  anode  substrate),  and  (v) 
an  RQ  element  (RiC:  diffusion  polarization  in  the  cathode),  which, 
however,  significantly  contributes  to  the  equivalent  circuit  only  at 
very  low  cathode  oxygen  partial  pressures  (p02  cat  <  0.05  atm),  cf. 
Table  1. 

At  this  point  it  is  needed  to  specify,  that  the  gas  diffusion  loss 
inside  the  7  |xm  thick  anode  functional  layer  (AFL)  can  be  neglected 
after  appositely  activating  the  cell:  The  AFL  is  a  dense  layer  after  ini¬ 
tial  reduction.  It  lasts  several  hours  (approx.  1 50  h)  of  cell  operation 
until  a  certain  open  porosity  in  the  AFL  structure  is  reached.  From 
there  on  gas  diffusion  in  the  AFL  can  be  neglected  for  cells  with 
an  AFL  thickness  smaller  than  8  p,m,  i.e.  for  the  cell  type  analyzed 
within  this  study  [7,9]. 

With  the  help  of  this  equivalent-circuit  model,  the  measured 
impedance  spectra  could  be  separated  into  the  single  loss  contribu¬ 
tions,  thus  permitting  a  precise  determination  of  the  area-specific 
resistance  values  for  the  single  loss  contributions  of  the  analyzed 
ASCs  (Table  1 ).  This  procedure  is  the  basis  for  the  experimental 
determination  of  the  modeling  parameters. 

4.  Results  and  discussion 

In  this  paper  the  I-U  characteristics  are  modeled  using  the 
Butler-Volmer  equation,  Fields  law  as  well  as  Ohm’s  law.  The  par¬ 
tial  pressure  dependency  of  the  exchange  current  densities  j0an  and 
jo, cat  of  anode  and  cathode,  respectively,  are  described  by  a  semi- 


Fig.  2.  Determination  of  parameters  c  and  d.  Ract,an  =R2a  +  R3a  as  a  function  of  pC 0 
(at  fixed  pC02  =  0.20  atm,  balance  N2 )  and  pC02  (at  fixed pCO  =  0.40  atm,  balance  N2 ) 
on  a  double-logarithmic  scale  (I  =  800  °C). 

empirical  power  law  ansatz.  The  temperature  dependency  of  the 
ohmic  resistance  and  of  the  exchange  current  densities  (jo, an  and 
jo, cat)  is  modeled  by  an  Arrhenius  ansatz.  The  complete  set  of  equa¬ 
tions  is  shown  in  Table  2.  For  more  details  on  the  model  the  reader 
is  referred  to  Ref.  [2]. 

In  order  to  make  a  prediction  of  the  anodic  activation  polariza¬ 
tion  overvoltage  (pact, an)  the  prefactor  yan ,  the  exponents  c  and  d, 
as  well  as  the  activation  energy  Eact,an  have  to  be  determined  (cf. 
equations  in  Table  2).  The  exponents  c  and  d  for  the  pCO  and  pC02 
dependency  are  determined  separately  in  this  study.  To  ensure 
that  pCO  and  pC02  could  be  varied  independently  of  each  other, 
a  defined  amount  of  inert  gas  (N2)  was  added  to  the  gas  mixture. 

In  Fig.  2  the  resistances  Ract,an(pC Oan)  and  Ract,an(pC02,an )  are 
plotted,  respectively,  as  a  function  of  partial  pressure  on  a  double- 
logarithmic  scale.  The  parameters  c= -0.058  and  d  =  0.25  can  be 
extracted  from  the  negative  value  of  the  slope  of  the  linear  fit,  in 
each  case.  The  negative  sign  is  a  result  of  the  inverse  proportionality 
between  exchange  current  density  and  the  resistance  (see  Table  2). 

To  determine  the  unknown  parameters  yan,  Eact,an,  B0hm 
and  Eact<ohm ,  impedance  measurements  were  carried  out  at  dif¬ 
ferent  temperatures  (650-800  °C).  The  ohmic  and  activation 
polarization  resistances  Rohm=R0  and  Ract,an=R2A  +  R3A  result¬ 
ing  from  the  CNLS  fit  are  shown  in  Fig.  3  as  a  function 
of  temperature.  These  resistances  were,  in  turn,  fitted  with 
the  Arrhenius  equation  (cf.  Table  2).  In  this  way  the  param¬ 
eters  yan  =  ( 4.56xl06T)Am-2,  Eact, an  =  118.64 kjmol-1  (1.23 eV), 


Table  1 

List  of  the  processes  known  to  take  place  in  a  ASC,  together  with  their  temperature  and  frequency  dependencies  as  well  as  the  characteristic  frequency  (Jr)  range  [2]. 
Process  fr  Dependencies  Physical  origin 

Pic  0.3-10  Hz  p02,  temperature  (low)  Gas  diffusion  in  the  cathode  structure  and  gas  flow  channel 

P2C  7-500  Hz  p02,  temperature  Oxygen  surface  exchange  kinetics  and  02--diffusivity  in  the 

bulk  of  the  cathode 

P-ia  0.070-10  Hz  pH2,pH20,pC0,  pC02,  temperature  (low)  Gas  diffusion  in  the  anode  substrate 

P2a  Psa  200 Hz-3  1<Hz3-50 kHz  pH2,  pH20,  pCO,  pC02,  temperature  (P2a+P3a)  gas  diffusion  coupled  with  charge  transfer  reaction 

and  ionic  transport  (AFL:  anode  functional  layer)  [1,7,9,11] 


Table  2 

List  of  model  equations  from  [2]. 
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Open  circuit  voltage  (OCV) 

Cell  voltage  (Uce„) 

Anode  and  cathode  gas  diffusion  overpotential 

{j]conc,el)  [13,14] 


Butler-Volmer  [15-19]  and  power  law  equations  [20] 
for  calculating  the  activation  overpotential  {r)actei) 

Relation  between  j0ei  and  Ract,ei 
Ohmic  overpotential  ( r}ohm )  [21,22] 
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Fig.  3.  Determination  of  parameters  yan,  Eact,an,  Rohm  and  Eacti0hm.  Arrhenius  plots  of 
Ract.an  =R2a  +  R3a  and  R0hm  =Ro- 

B0hm  =  5.8  x  1012  SKnrr2  and  Eactohm  =  91.66  kj  mol-1  (0.95  eV)  are 
obtained  (cf.  Fig.  3  and  Table  3). 

The  structural  parameters  xj/an  =  0.133  and  xl/cat  =  0.022  for 
the  calculation  of  the  gas  diffusion  polarization  and  the 
parameters  ycat  =  (1-52  x  108T) Am-2,  m  =  0.22,  acat  =  0.65  and 
Fact, cat =  139.86  kj  mol-1  for  the  calculation  of  the  cathodic 
exchange  current  density  j0<cat  were  already  determined  in  a  pre¬ 
vious  study  (on  nominally  identical  cells),  and  thus  adopted  here 
[2].  Therefore,  the  sole  remaining  parameter  to  be  determined 
is  the  charge  transfer  coefficient  aan •  Different  to  the  previous 
procedure,  this  parameter  is  obtained  by  impedance  measure¬ 
ments  under  load,  since  the  charge  transfer  coefficient  affects  the 


Temperature,  T/°C 
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Table  3 

Modeling  parameters  determined  from  impedance  measurements. 


Parameter 

Value/expression 

Unit 

Comment 

Rohm 

5.8  x  1012 

SKirr2 

Eact,ohm 

91.66 

kj  mol-1 

C 

-0.058 

- 

d 

0.25 

- 

m 

0.22 

- 

Adopted  from  [2] 

Van 

4.56x  10 6  xl 

Am-2 

Ycat 

1.52  x  10s  xf 

Am-2 

Adopted  from  [2] 

Eact,an 

118.64 

kj  mol-1 

Eact,cat 

139.86 

kj  mol-1 

Adopted  from  [2] 

Oian 

0.62 

- 

Meat 

0.65 

- 

Adopted  from  [2] 

an 

0.133 

- 

Adopted  from  [2] 

cat 

0.022 

- 

Adopted  from  [2] 

current-dependent  behavior  of  the  activation  overpotential.  aan 
was  obtained  by  a  (polynomial)  fit  of  the  evaluated  current  depen¬ 
dent  Ract,an  followed  by  a  least-squares-fit  of  the  Butler-Volmer 
equation  (for  details  on  this  procedure  the  reader  is  referred  to 
Ref.  [2]).  The  value  determined  in  this  way  {aan  =  0.62)  is  listed  in 
Table  3. 

To  validate  the  model,  I-U  characteristics  were  recorded  at  dif¬ 
ferent  temperatures  and  gas  compositions.  In  Fig.  4a  and  b  the 
measured  and  predicted  I-U  curves  are  compared.  As  can  be  seen, 
the  model  delivers  a  very  reliable  prediction  of  the  measured  curves 
over  the  entire  parameter  range. 

At  high  operating  temperatures  the  electrochemical  reactions 
at  the  electrodes  are  faster,  as  a  result  the  overpotentials  due 
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Fig.  4.  Simulated  (lines)  and  measured  (symbols)  I-U  curves  at  four  different  tem¬ 
peratures  (a)  and  gas  compositions  (b). 
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to  activation  polarization  are  lower.  Moreover,  the  temperature 
dependence  of  the  ohmic  losses  R0hm(T )  exhibit  an  Arrhenius-type 
behavior,  which  also  leads  to  lower  ohmic  losses  at  high  tempera¬ 
tures.  However,  higher  temperatures  also  yield  lower  open-circuit 
voltages.  These  temperature-dependent  characteristics  typical  of 
the  I-U  curves  can  be  observed  for  all  predicted  curves  (Fig.  4a). 

By  analyzing  Fig.  4a  it  can  be  seen  that  the  model  traces  beyond 
the  cell  voltage  of  0.6  V  are  approaching  a  limiting  current  den¬ 
sity  (at  Uceu  =  0)  which  is  weakly  dependent  on  temperature.  This 
behavior  can  be  explained  by  the  fact  that  the  gas  diffusion  polar¬ 
ization  is  characterized  by  a  negligible  dependency  on  temperature 
as  already  discussed  in  Refs.  [1,2,13]. 

Fig.  4b  shows  the  simulation  results  for  different  CO  and  CO2 
partial  pressures  at  a  constant  temperature  (T=800°C)  and  oxy¬ 
gen  partial  pressure  at  the  cathode  (pO2)Cat  =  0.21  atm  (air)).  Low 
CO  partial  pressures  were  approached  appositely  to  validate  the 
model  within  the  diffusion  limitation  regime  of  the  I-U  charac¬ 
teristics.  For  the  cell  type  investigated,  gas  diffusion  polarization 
at  the  anode  becomes  already  predominant  at  carbon  monoxide 
partial  pressures  pCOan  <  0.4  atm.  A  comparison  between  simulated 
and  measured  I-U  curves  at  pCOan  =  0.24  atm  clearly  shows  that 
the  diffusion  overpotential  at  the  anode  predicted  by  the  model  is 
sufficiently  precise. 

5.  Conclusions 

A  zero-dimensional  stationary  model  for  the  prediction  of 
current-voltage  characteristics  of  planar  anode-supported  SOFC 
single  cells  operated  on  CO  has  been  presented.  Simulation  results 
show  a  very  good  agreement  with  experimentally  determined 
data.  The  model  is  based  on  the  Butler-Volmer  equation  for  the 
determination  of  activation  polarization,  Fields  model,  using  the 
calculation  of  the  structural  parameters  ( ij/an  and  i/fcat),  for  obtain¬ 
ing  the  diffusion  polarization,  and,  lastly,  Ohm’s  law  for  calculating 
the  ohmic  overpotential.  The  necessary  modeling  parameters  were 
all  obtained  from  impedance  measurements,  thereby  specifically 
varying  the  operating  conditions,  e.g.  fuel  composition  at  the  anode, 
oxidizing  gas  composition  at  the  cathode,  current  load  density  and 
temperature. 

However,  the  model  only  yields  precise  results  assuming  a 
spatially  constant  current  density  and  temperature  distribution 
throughout  the  active  electrode  area  which  requires  cells  with 


a  sufficiently  small  cross-sectional  area  (1  cm2  electrodes  in  our 
case).  Moreover,  in  the  future  the  model  shall  be  employed  for 
the  prediction  of  the  performance  of  larger  single  cells  (stack- 
“repeating  units”).  For  this  purpose,  it  has  to  be  expanded  into  a 
two-dimensional  model  in  order  to  consider  a  spatially  resolved 
current  density  and  the  resulting  in-plane  temperature  distribu¬ 
tion  of  the  cell.  Of  course,  in  case  of  larger  cells  the  voltage  drop 
due  to  the  fuel  utilization,  which  has  been  neglected  so  far,  has  to 
be  considered. 
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